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We analyze non-equilibrium current transport in molecular electronic devices, using as an example 
devices formed by two terphenyl dithiol molecules attached to gold electrodes. Using a first-principles 
based self-consistent matrix Green's function method, we show that the spatially resolved current 
density profiles and local electrochemical potential drops provide valuable information regarding the 
local field effect on molecular transport, which depend on the internal structure of the molecules 
and cannot be obtained from measuring the current- and conductance- volt age characteristics alone. 
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Introduction. — Detailed understanding of electron 
transport through single molecules contacted by metallic 
electrodes is crucial for the development of the emerg- 
ing technology of molecular electronics. 1 Although much 
attention has been devoted to the measured molecular 
conductance and its dependence on metal-molecule in- 
teractions, 1 less is known about the non-equilibrium as- 
pects of the phase-coherent molecular transport, which 
provide information not obtainable from the conductance 
measurement. 234 

Several important effects are of interest for non- 
equilibrium transport: (1) A finite electric field asso- 
ciated with the applied bias voltage will induce charge 
pileup (screening) within the molecular junction. The lo- 
cal transport field experienced by the tunneling electron 
is the self-consistent screened field; (2) Since the total 
current is conserved throughout the molecular junction 
for dc-transport, the spatial distribution of the current 
density can be highly non-uniform since the molecules 
are intrinsically inhomogeneous; (3) As current flows, the 
electrons at the source and drain contacts have differ- 
ent electrochemical potentials. It is not clear if an effec- 
tive local electrochemical potential (LEP) can be defined 
everywhere within the molecular junction to character- 
ize the electron distribution in the non-linear transport 
regime. 

The screened nature of the local transport field has 
been analyzed in previous worksj2i4 following early ideas 
of Landauer and Buttiker in mesoscopic transport sys- 
tems^ In particular, resistivity dipoles and strong local 
fields have been found in the vicinity of potential bar- 
riers both at the interface and at certain points inside 
the molecules. These in general modify both the en- 
ergy and wavefunction of the frontier molecular states 
and lead to different electrostatic potential drops depend- 
ing on the molecules and the device structures.^ The 
non-equilibrium electron distribution within the molec- 
ular junction has also been analyzed by projecting the 
nonequilibrium density-matrix into individual molecu- 
lar orbitals^ Here, we extend previous work by present- 
ing real-space analysis of current density profiles and lo- 
cal non-equilibrium electron distributions throughout the 



molecular junctions. Such analysis has important impli- 
cations in resolving the following issues: (1) The identi- 
fication of current path gives a clear picture of the con- 
ducting and non-conducting part of the molecular junc- 
tion. In particular, large current density may develop 
at certain parts of the molecule, which forms a bottle- 
neck to transport and might possibly be imaged using 
scanned probe microscopes; 6 (2) The local electrochemi- 
cal potential (LEP) characterizes the energy-distribution 
of the electrons everywhere within the junction, and is 
relevant in multi-probe transport measurement using a 
weakly coupled phase-sensitive voltage probe, 7 e.g., in 
scanning tunneling potentiometry (STP). 8 The LEP de- 
fined here is determined by the local density of states 
(LDOS) of the molecular junction in the presence of cur- 
rent transport, which may possibly be probed using scan- 
ning tunneling spectroscopy (STS)AiS 

N on- Equilibrium Green's Function (NEGF) approach 
to scattering theory of molecular transport. — We ana- 
lyze the non-equilibrium effects within the scattering the- 
ory of mesoscopic transport^ where the relevant phys- 
ical quantities are obtained from microscopic calcula- 
tions using a self-consistent matrix Green's function the- 
oryjiiS which combines the NEGF theory™ with an ef- 
fective single-particle description of molecular junction 
electronic structure using density functional theory^* We 
use atomic unit throughout the paper unless otherwise 
noted. 

The central quantity within the NEGF theory 
is the correlation Green's function G < (f,r';E) = 
'^ijGfj{E)4 l i{r)<l)j(r') 1 where we have expanded G < 
in terms of local atomic basis sets <pi (which are 
real functions)^ The current density is obtained from 
G < (f,r';E) asi£ 

J(r) = 1/2 limJV'-V) / dEG < (f,f':E) =l/2^( / dEG 

i,j 

(1) 

where we used the notation dSij(r) = <j>i{r)V(j)j{r) — 
(f)j(r)\J 4>i(r). The terminal current is obtained 
by integrating the current density over a sur- 
face separating the molecule from the con- 
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FIG. 1: (Color online) (a) Schematic illustration of the 
metal-molecule-metal junction for the two terphenyl dithio- 
late (TPD) molecules. Six gold atoms on each metal sur- 
face (12 overall) are included into the "extended molecule" 
where the self-consistent calculation is performed. The effect 
of the rest of the electrodes (with the 6 atoms on each side 
removed) are modeled as self-energy operators. Also shown is 
the bias polarity and coordinate system of the molecular junc- 
tion; (b) The left figure shows the equilibrium electron trans- 
mission coefhcent. The right figures show the self-consistent 
current- voltage (I-V) and differential conductance- voltage (G- 
V) characteristics of the two molecular junctions. 



tact, which leads to the standard formulaii*i2ii£ 

Il { b) =JdETr{T HR) [f(E- t x HR) )A(E)+iG<(E)]} = 
JdEjdf{T L{R) (f,r)j: tj [f(E fi L{R) )A(f,r;E) + 
iG<(r,r;E)]}, where A(r,r;E) = ^ A^E^rj^r). 
The spectral function A and the operators ^l(r) de- 
scribing the contacts are defined in the standard 
manner Note that in the matrix notation, 
we have A = G R [T L + T R ]G A = A L + A R and 
~iG< = G R [T L f L + T R f R ]G A = A L f L + A R f R , 
where A L ^ R ) = G R T L ^G A . Defining a lo- 
cal transmission function for electron injection 
from the left (right) electrode TumffiE) = 



FIG. 2: (Color online) Cross sectional view of current density 
profiles at gold- TPD 1-gold (a) and gold-TPD2-gold (b) junc- 
tions. Also shown is the position of the molecules. For each 
junction, the upper figures show the direction of current den- 
sity, the lower figures show the magnitude of current density 
component along the transport direction (V axis here). 



T L ( R )(f,r;E)A(r,r;E), the total current can be 
expressed in terms of T L ^ R ^(f; E) and an effective 
local electron-distribution function f e ff(f;E) as: 
Il ( r) = fJdEj drT L(R) (r; E) [f L{R) (E) - f eff (r; E)] , 
where A(r,r; E)f ef f(r; E) = -iG < (r,r:E) = 
A L (r, r; E)f L {E) + A R (r, r; E)f R {E). 

Replacing the spectral functions by the corresponding 
density of states A(r, r; E) = 2np(f; E), ^l(a) r; E) = 
27r/9£(^j)(r; E), we arrive at the following definition of an 
effective local non-equilibrium electron distribution func- 
tion: f eff (r;E) = ^§.f L {E) + £^f R (E). At low 
temperature and within the linear-response regime, we 
can replace the Fermi distributions by step- like functions 
and define the local electrochemical potential (LEP) as 

Me//(f) = E W§jy Mi + P p^eI) ^R where E f is the equi- 
librium Fermi-level. At finite temperature in the non- 
linear transport regime, the above definition of LEP re- 
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FIG. 3: (Color online) Cross sectional view of local elec- 
trostatic potential and electrochemical potential drops at the 
gold-TPDl-gold (a) and gold-TPD2-gold (b) junctions. Also 
shown are the positions of the molecules. The electrostatic 
potential varies rather smoothly from —V/2 to V/2 across the 
junction, but the electrochemical potential oscillates strongly 
between -V/2 and V/2. 



mains useful for understanding transport physics if the 



coefficients 
[iL and jj,R 



PL{R)(r;E) 



vary slowly with energy between 



p(r;E) 

Since this is true for molecules chcmisorbcd 



onto the metal surfaces, we get the following LEP defined 
everywhere within the molecular junction: 



measured by an ideal phase-sensitive non-invasisve volt- 

In 



Veff(r) 



Ejj PL;ij{Ef)(j>i{r)<j)j{r ) i T,ijPR;ij{ E f) ( l>i( r l < P. 
Ey Pij( E f)<i>i(r)<l>j(rj 



"Mi" 



(2) 

Note that an equivalent definition of the local electro- 
chemical potential has been given by Gramespacher and 
Biittikcr 15 for a mesoscopic conductor within the scat- 
tering matrix theory. Here the contact-resolved local 
density-of-states Pl(r) = G R T l(r)G a /27r is proportional 
to the injectivity of the left (right) electrode given by 
Gramespacher and Biittiker. The LEP /x e //(r) as defined 
here in a current-carrying conductor may possibly be 



age probe (Biittiker probe) introduced at point filH 
scanning tunclling potentiometry, this is identified with 
the bias voltage between the STP tip and the sample 
under the condition of zero tunneling currenti^ 

Results and their interpretation. — Since local-field ef- 
fects are most clearly seen for longer molecules with 
internal barriers^ we apply our theory to the devices 
formed by two extended tt systcms-a three-ring oligomer 
of phenylene ethynylene with dithiol substituents (called 
TPDI) and a terphenyl dithiol molecule (called TPD2)- 
in contact with gold electrodes through the end sul- 
fur atoms. The device structure is illustrated schemat- 
ically in Fig. ma). The benzene rings are co- planar 
when connected by a triple-bonded C-C bridge (TPDI), 
which leads to optimal orbital overlap. In the absence 
of the C-C bridge (TPD2), a torsion angle of 36° is in- 
duced between neighbor benzene rings, which leads to 
weaker orbital overlap and an effective potential barrier 
for electron motion inside the molecule £ The calcula- 
tion is performed using a modified version of Gaussian98 
program^ with the Becke-Perdew-Wang parameteriza- 
tion of density-functional theorjiii and appropriate pseu- 
dopotentials with corresponding optimized Gaussian ba- 
sis sets. The details of the self-consistent calculation have 
been discussed extensively elsewhere All 

The equilibrium (zero bias) electron transmission char- 
acteristics of the molecular junctions are shown in Fig. 
I2b), from which we find that for both molecules, 
the highest-oocupied-molecular-orbital (HOMO) lines up 
closer to the metal Fermi-level Ep than the lowest- 
unoccupied- molecular-orbital (LUMO), similar to previ- 
ous finding on phenyl and biphenyl dithiolate molecules, 3 
The lineup scheme is more favorable for the pla- 
nar TPDI molecule. The room-temperature con- 
ductance of the molecular junction is 0.82(/iS f ) and 
OM(nS) for TPDI and TPD2 junctions respectively. 
The self-consistent current-voltage (I-V) and differen- 
tial conductance- voltage (G-V) characteristics in the bias 
range of — 2(V) to 2(F) are shown in Fig. IHb). Note that 
despite the > 0.3(eF) difference in the HOMO level po- 
sition with respect to Ep at equilibrium, there is only a 
0.1(V) difference in the bias voltage where the junction 
reaches the first conductance peak (1.2(F) and 1.3(V) 
for the TPDI and TPD2 junction respectively). This is 
due to the stronger bias-induced modification of molec- 
ular states (static Stark effect) in the TPD2 molecule, 
e the internal barrier across the neighbor benzene 
stronger local-field variations inside the 
This pushes the HOMO up relative to 



;s fl leads 



to 

molecule core, 
the equilibrium Fermi-level and correspondingly reduces 
the voltage needed to move the metal Fermi-levels past 
the HOMO level. 

The overall I-V and G-V characteristics of the two 
molecules show similar behavior, but the spatially- 
resolved current density distributions are quite differ- 
ent due to the different molecular structures, which are 
shown in Fig. [21 for the two molecules at bias voltage 
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of 2(V). Note that the direction of current is reverse 
to that of the net electron density flux due to the neg- 
ative electron charge. For the planar TPD1 molecule, 
we show both the vectorial plot (direction) of current 
density and the magnitude of the y-component (perpen- 
dicular to the electrode surface) of current density j y in 
the xy-plane (defined by the benzene rings) and in the 
plane located 1.5(A) above. For the nonplanar TPD2 
molecule, this is shown in the xy-plane (defined by the 
left-most and right-most benzene rings) and in the plane 
located 2.5(A) above (the carbon and hydrogen atoms in 
the central benzene ring are located at Z = ±0.7 (A) and 
±1.3 (A) respectively). 

For the planar TPD1 molecule, the small y-component 
of the current density in the xy plane peaks around 
the peripheral hydrogen atoms of the central benzene 
ring. j y increases rapidly (and symmetrically) moving 
above/under the benzene plane, peaks near 0.8(A) from 
it and then decreases slowly with further increase of \Z\. 
Note that the peak location of j y corresponds to the peak 
location of the density of the it electrons in the molecule. 
The direction of the dominant (7r-electron) current flow is 
such that electrons are injected from the right electrode 
mostly along the sulfur-surface bond, and then propa- 
gate through the middle of the right-most benzene ring 
(centered around Y"-axis), being diverted to the perime- 
ters of the central benzene, and converge again into the 
middle of the left-most benzene ring. For the nonplanar 
TPD2 molecule, there is a large j y in the plane Z = 0, 
which resides mostly around the central benzene ring. 
Moving away from Z = 0, j y increases gradually in both 
the left and right benzene rings, and reaches its peak at 
about 0.7(A), corresponding to the position of the car- 
bon atoms in the central benzene (the peak value is only 
about three times that in the XY plane) and then de- 
creases slowly. For positive (negative) Z, j y is larger on 
the negative (positive) X half of the junction, determined 
by the location of carbon and hydrogen atoms in the cen- 
tral benzene ring. The direction of the dominant current 
flow is such that electrons propagate rather uniformly 
through the right-most benzene ring, being diverted to- 
wards the negative (positive) X half of the molecule at 
positive (negative) Z and propagate again rather uni- 
formly through the left benzene ring. 

The local electrostatic potential and electrochemical 
(LEP) potential drops in the molecular junction are 
shown in Fig. [3] for bias voltage of 2(V). Unlike the 



electrostatic potential drop which varies rather smoothly 
from — V/2 to V/2 across the junction, the local elec- 
trochemical potential shows oscillatory behavior between 
—V/2 and V/2 throughout the molecular junction. This 
is due to the phase-sensitive nature of the voltage probe: 
an electron wave incident from either electrode gives two 
contributions to the electron flux in the voltage probe at 
point f, i.e., the direct transmission and the transmis- 
sion after multiple scattering within the molecular junc- 
tion. The oscillatory behavior in LEP therefore demon- 
strates the presence of a potential barrier for electron 
injection into the molecular junction. The more opaque 
(smaller transmission coefficient) the barrier is, the larger 
the magnitude of the oscillation will bei The oscillation 
of the LEP is also sensitive to the shape of the poten- 
tial barriers, which shows up clearly in the two planes of 
the molecular junctions (Fig. 0. The LEP in the TPD2 
junction oscillates more strongly than the TPD1 junction 
due to the presence of barriers both at the metal- molecule 
interface and inside the molecule induced by the weaker 
orbital overlap across neighbor benzene rings. In addi- 
tion, there is also oscillation in the direction (along X 
axis) perpendicular to the transport direction due to the 
three-dimensional structure of the molecular junction. 

Discussion and Conclusion. Although in principle 
both the local electrochemical potential (LEP) and the 
current density profile in single-molecule devices may 
possibly be measured using scanning nanoprobe tech- 
niques, in practice this can be extremely difficult due to 
the spatial resolution involved and the requirement that 
the local probe should be minimally invasive. For exam- 
ple, the LEP varies on the scale of the Fermi wavelength 
(\f) of electrons injected from the electrodes, which is ~ 
1(A) for metallic electrodes, unlike the two-dimensional- 
electron-gas fabricated from semiconductor heterostruc- 
tures where Xp can be of 100(nm) and longer^ Fur- 
ther work is therefore needed for a quantitative evalua- 
tion taking into account realistic experimental conditions 
in such multi-probe transport measurements, in partic- 
ular the probe geometry effect and local probe-induced 
perturbation of the junction electronic processes. The 
present work hightlights the importance of local atomic- 
scale analysis in revealing subtle effects in single-molecule 
electronic devices. 
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